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ABSTRACT 

The nearby Type la supernova SN 2011fe in M101 (cz = 241 km s _1 ) provides a unique opportunity 
to study the early evolution of a "normal" Type la supernova, its compositional structure, and its 
elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first 
month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a 
clear picture of how various line-forming species are distributed within the outer layers of the ejecta, 
including that of unburned material (C+O). We follow the evolution of C II absorption features until 
they diminish near maximum light, showing overlapping regions of burned and unburned material 
between ejection velocities of 10,000 and 16,000 km s _1 . This supports the notion that incomplete 
burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among 
SNe la. The observed evolution of the highly Doppler-shifted O I A7774 absorption features detected 
within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 
21,000 km s~F The fact that some O I is present above C II suggests that SN 2011fe may have had 
an appreciable amount of unburned oxygen within the outer layers of the ejecta. 

Subject headings: supernovae: general supernovae: individual (SN 201 lfc) 
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1. INTRODUCTION 

SN 2011fe in M101 is the closest (6.9 Mpc) Type la 
supernova (SN la) in 25 years. It was discovered only 11 
hours after explosion, the earliest of any SN la (Nugent 



et al. 2011) and has revealed more about its progenitor 
than any other recent SN la. In the canonical picture 
of a SN la, a CO white dwarf (WD) gains mass until 
it approaches the Chandrasekhar mass, either by mass 
transfer from a nondegenerate comp anion, or by m erging 
with another WD (for a review, see |Howell||20lTl ). 

The merger of two WDs has been strongly suspected 
in the case of "super-Chandra" SNe la, those that are so 
luminous they appear to require a progen i tor above the 
Chandrasekhar mass (Howell et al. 2006 Scalzo et al. 



20101. Unlike most SN la, these subtypes are known tor 
showing conspicuous and relatively long-lasting spectro- 
scopic signatures of carbon in their outer la yers, possi- 
bly left ove r from the violent CO- WD merger ( Silverman 
et al.|[20lT |. But this scenario cannot explain all SNe la 



as some show narrow lines of sodium from circumstellar 
material (CSM), thought to be from a nonde generate sec- 
ondary star in a single- degenerate system (|Di lday||2012 
and references therein). [Sternberg et al. ( |201ip estimate 
that at least 20% of SNe la in spiral galaxies arise from 
this channel. Furthermore, it remains to be seen if cir- 
cumstellar sodium is concurrent with carbon detections. 
Could abundance studies of carbon lines reveal the pro- 
genitor system? 

With SN 201 lfe, it was confirmed observationally for 
the first time tha t the exploding star was a WD ( |Nu-| 
gent et al. 2011 Bloom et al. 2012). Pre-explosion 
Hubble Space Telescope images revealed that any sec- 
ondar y star could not be a luminous red giant (Li et al. 
20111. Furthermore, the lack of an observed shock, ex- 



pected when the ejecta hits any nondegenerate secondary 
(Kasen 2010 1, suggests that the secondary could not be 
a red giant ( Nugent et al. 2011 ), and p robably not any 
nondegenerate star (Bloom et al. 2012). Nondetections 
of radio and X-ray emission place strong limits on the 
amount of CSM, disfavoring most single-degenerate hy- 
potheses, but consistent with the merger of two WD stars 



([Horesh et al.|2012|[Chomiuk et al.|2012||Margutti et al 
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In addition, thorough and complete spectroscopic ob- 
servations of SN 201 lfe offer a chance to study the dis- 
tribution of ejected material. For instance, as many as 
a third of SNe la show evidence of unburncd material if 



spectra are obtained early enou 


gh (Parrent et al. 2011 


Thomas et al. 


|2011a| IFolatclli etahl 


2012; Silverman k 


Filippcnko||2012| |Blondin et al.||2012 


). Carbon features 



( Scalzo et al.||2010" ), but rather signal some difference in 
the e xplosion physics of th e primary WD (Baron et al 



2003| |Thomas et al.||2007[ ). Abundance studies of un 
burned material can be coupled with constraints on ma- 
terial that is synthesized in the explosion, as some ex- 
plosion models predict the co mplete burning of the WD, 
leaving no u nburned carbon (Kasen et al. 2009 Maeda 
eFlLl|2010b ). 
To invest ie 



gate these questions it is necessary to com- 
bine our limited progenitor knowledge with a map of the 
elemental distribution of ejected material in SNe, includ- 
ing unburned material. Here we present 18 high signal- 



to-noise (S/N) optical spectra of SN 201 lfe during its 
first month, starting one day after explosion. We fit 
all spectra with the automated spectrum synthesis code, 
SYNAPPSp 4 } to trace the compositional structure of the 
ejected material. 

2. OBSERVATIONS 

On 2011 August 24.167 (UT dates are used through- 
out this Letter), the Palomar Transient Factory (PTF) 
discovered SN 201 lfe in M101. The SN was discovered 
early enough that the time of explosion, t cxpi could be 
estimated to within ~ 20 minutes of 2011 August 23.69 
The SN reached a -B-band pea k 



(Nugent et al. 20111 



brightness on September 10.3 ( Bianco||in preparation 
Optical spectra were obtained on eight different tele 
scopes (Table 1, Figure [I]). Reductions were done with 
IRAF23 

3. SPECTROSCOPIC ANALYSIS 

From the onset, the spectroscopic features of SN 
201 lfe appear broad and comprised of P-Cygni ab- 
sorption/emission profiles of intermediate-mass elements 
(IMEs) and iron-p eak elements (IP Es) typically seen in 
super nova spectra ( Branch fc Patchett||1973| |Filippenko 
1997[ ). Over time the features narrow, depicting an evolu- 
tion similar to "normal" subtype objects. Most notably 
the spectra exhibit: (1) high velocity features (HVFs; 
v > 16,000 km s^ 1 ) of Si II, Ca II and Fe II absorption 
lines from the onset of the explosion, (2) unburned ma- 
terial in the form of a C II A6578 absorption feature, and 
(3) O I A7774 at velocities higher th an previously de- 
tected for a SN la ( |Nugent et al.|2011| . From the lack of 
conspicuous Ti 11/ be III absorption features sometimes 
observed in sub- and over-luminous SNe la, respectively, 
as well as the slowly decreasing blueshift of the Si II 6150 
A feature, our spectroscopic observations of SN 201 lfe 
in Figure [I] are consistent with characteristics of a "low 
velocity gr adient" / "core normal" subtype (see Blondin 
|et al.|[2012] and references therein). 

3.1. Line Identifications: SYNAPPS 

Given that supernova expansion velocities are ^10 4 km 
s , Doppler line blending can be an obstacle for map- 
ping the underlying compositional structure. SYNAPPS 
is a useful tool for identifying line features in photo- 
spheric phase spectra through its empirical spectrum 
synthesis procedure which requires no a priori composi- 
tional structure as input. Under the assumption of pure- 
resonance line transfer and a sharp blackbody-emitting 
photosphere, an atomic line list is used to construct 
a basic spectrum from a pseudo-blackbody continuum 
level. The optical depth profile for a line is modeled as 
t(v) = r re f x exp(^ r ^ — where v is the expansion ve- 
locity, v re f is an arbitrary normalization velocity, v e is 
the e-f olding velocity, and r re / = r(v re f) (Thomas et al. 
2011b I. 



24 Derived from SYN0W, see iBranch et all <|2005|) and iThomasI 
|et^j2UlTEt . 

2U iRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (NSF). 
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Because the sources of line formation influence nearly 
all regions of the spectrum from a wide range of depths 
within the ejecta, we model the spectrum approximately 
by treating the minimum reference velocity u m i n , for ev- 
ery ion as a free parameter. For a given set of ions in a fit, 
we set the photospheric velocity w p hot > equal to the low- 
est measured v m i n . This follows the notion that not all 
blended absorption minima coincide with the Doppler 



tively infer the presence of ions such as cobalt and nickel 
from the first month of optical observations alone, i.e. 



shifted loca tion of the photosphere over time (Bongard 
et al. 2008 ) . Furthermore, we use a quadratic warp- 
ing function to simplify the treatment of the continuum, 
calibration errors, and uncorrected reddening. Conse- 
quently, our fits are constrained by how we are able to 
best match absorption minima and the relative strengths 
of all features. 

For the spectra in Figure [T[ we converged a series of 
SYNAPPS fits. In Figure[2]we plot the spectra for the days 
1.5, 10.4, 20.6 and 32.6 post-explosion, as they best doc- 
ument the relevant temporal changes on which we focus 
our attention. For each epoch, we compare our SYNAPPS 
fit and over-plot the continuum fit for reference. We have 
loosely highlighted various regions of the fits to indicate 
which parts of the data we use to infer, and thereby con- 
strain, the presence of particular ions. For example, the 
two broad absorption features centered around 4300 and 
4800 A are complex blends of three or more ions. When 
we remove any one of these ions from the fit, we are un- 
able to compensate for it with any reasonable perturba- 
tions of the remaining parameters. That is to say, while 
one or two ions may dominate it, the shape of a par- 
ticular blended feature is still subject to the minor but 
observable influences of other species. 

For the 1.5 day fit we set v p hot = 16,000 km s" 1 and in- 
clude photospheric components of C II, O I, Mg II, Si II, 
S II, Ca II, and Fe III. In addition, we include higher 
velocity components of O I, Si II, Ca II, and Fe II above 
16,000 km s _1 in order to follow the evolution of the cor- 
respondi ng HVFs detected during the earliest observed 
epochs ( |Nugent et al. 2011 1. Overall, many of the fea- 
tures are accounted tor m our subsequent fits, and the 
relative fluxes and respective minima match nicely. 

Over the course of their detection, we model the Ca II 
HVFs as material "detached" 10,000 km s _1 above the 
photospheric line forming regions. It is clear when the 
Ca II HVFs disappear for SN 201 lfe at ~5 days after 
texpi] however, the same is not true for Si II. We model 
Si II as present in both photospheric and higher velocity 
regions, despite the separation of 2000 km s _1 near max- 
imum light. It is sometimes clear that the 6150 A fea- 
ture requires two distinctly separate compo nents of Si II 



for a best fit (most notab ly SN 2009ig ; see Foley et al. 
2012|. Recently however, Smith et al. (20111) presented 



spectropolarimetry of SN 201 lfe showing two Si II ab- 
sorption line polarizations at ^7300 and 12,000 km s _1 
just after maximum light (see their Figure 3). During 
this phase our two components of Si II are at 8000 and 
11,000 km s _1 . While the overlap is not exact, the find- 
ing of a two-component Si II feature is consistent with 
the spectropolarimetry. 

The set of ions mentioned above is primarily the same 
used until maximum light, at which point Na I, S II, and 
Fe II begin to account best for most of the features. For 
SN 201 lfe, while using only SYNAPPS, we cannot dcfini- 



species responsible fo r powering SN la lightcurves (Col- 
gate & McKee 1969). However, early-phase broadband 
photometry, as well as spectros copic analyses of infrared 
data, could help in this regard ( Marion et al.||2009 Piro 
20121. 



3.2. Spectroscopic Modeling: SYNAPPS 

The result of any series of SYNAPPS fits is the esti- 
mation of velocities spanned by a particular set of ions. 
Therefore, we can use spectra in two different ways, as 
seen in Figure|3] From our time-series fits, in Figure|3ja) 
we show our v m i n (t b s ) estimates for each of the pro- 
posed ions in §3.1. A clear separation between outer 
and inner regions of ejected material can be seen. The 
material near photospheric velocities that span roughly 
7000 - 16,000 km s" 1 includes Na I, C II, O I, Mg II, 
Si II, S II, Ca II, and Fe III, whereas most of the higher 
velocity material between 16,000 and 30,000 km s" 1 is 
comprised of O I, Ca II, Si II, and Fe II. 

Since we treat u m m as a free parameter, there is some 
dispersion among ions near the modeled photosphere for 
a given epoch (black dashed line in Figure [3l a)). How- 
ever, SYNAPPS fits are generally good to within 250 km 
s _1 and the scatter of our inferred v m i n is no less than 
500 km s . Therefore we consider the dispersion con- 
sistent with what should be expected in an environment 
of multiple line forming regions. 

When we take into account the Doppler-widths of our 
fits, we are able to place tighter constraints on the lo- 
cation of species by constructing a velocity-scaled map 
that traces how far each ion extends. In Figure [3jb) we 
plot a single normalized absorption profile for each ion, 
summed over the first month. Clearly, a photospheric 
and higher velocity region can be seen. In addition, Fig- 
ure j3^b) yields some insight into where species are likely 
to be present, with respect to one another, throughout 
the ejecta. This is particularly important for investigat- 
ing unburned material since its location with respect to 
burned material can be used to constrain the details of 
the explosion nucleosynthesis. 

Figure [4] highlights regions of interest for tracing un- 
burned material with C II A6578 and potentially O I 
A7774 (see also Figure 3]). For our C II fits, we are able 
to follow the suspecteaabsorption feature as it evolves 
across the emission component of the Si II A6355 line 
over time. This implies that unburned carbon is present 
at least between 10,000 and 16,000 km s _1 . Our O I fits 
in Figure [4] require both photospheric and higher velocity 
components of O I for the 1.2, 1.5, 2.6 and 4.5 day spec- 
tra as the O I HVF disappears quickly and within 5 days 
after t exp . Within the time-frame of our observations, 
C II A6578 is present from 1.5 to 20.6 days after t exp , 
i.e., from the onset and until at least maximum light. 
The same is true for O I; however we find that O I ini- 
tially extends at least 5000 km s _1 above C II, within the 
outermost layers. This suggests the O I detected during 
the earliest epochs is likely to be unburned progenitor 
material rather than a product of carbon burning in the 
SN. 

4. DISCUSSION & CONCLUSIONS 
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Thanks to the proximity of SN 2011fe, we have 
been able to obtain frequent high S/N spectra and use 
SYNAPPS to produce what may be the best-ever map of 
the material distribution in a "normal" SN la for the first 
month of its evolution. For SN 201 lfe, it is important 
to note that, apart from differences in relative expansion 
velocities (Figure [3]ja)), the spectra do not vary signif- 
icantly from one epoch to another. Furthermore, as in 
most SNe la, in SN 201 lfe we find O I, Mg II, Si II, S II, 
Ca II, and Fe II present from 16,000 km s _1 or more 
down to about 7000 km s _1 (later spectra will probe 
deeper). This implies similar compositional constituents 
over a large region of ejected material and is supported 
by our SYNAPPS fitting results (Figure [^b)). 

Stretching to nearly 30,000 km s , high velocity 
regions of silicon and calcium are inferred from pre- 
maximum light spectra. A region of sodium is inferred 
to be present at high velocities as well, but is detected 
post-maximum light when the outer layers begin to cool. 
However, rather than l arge-sc ale asymmetries, as pre- 
dicted by Maeda et al. (2010a), these high velocity fea- 
tures may trace only a small a mount of material, as little 
as 10~ 7 Mf7) in t he case of Ca II (j H atano et al.|l999||Maz-| 
zali et al. 2005 Tanaka et al. 2005). Such a small inferred 
Ca mass supports the notion that strong features in SN 
spectra do not always correspond to large masses of line- 
forming ejecta. 

We also find Si III and Fe III present at early times (be- 
fore 3 days after maximum light), indicating that some 
parts of the ejecta achieved high temperatures, although 
the entire photosphere was not excessively hot. Even 
a day after explosion the spectra were not particularly 
blue, nor were they dominated by high-ionization species 
Various scenarios that predict an early shock fr om in- 



terac tion with a nondegenerate secondary star (Kasen 



2010|, a torus of m aterial from a disrupted secondary 
(|Scalzo et al. 20101, or CSM, are not supported by the 



data. This is consistent with a lack of radio a nd X 
ray detections ffloresh et al.||2012 Chomiuk et al 



2012 



Margutti et al.||2"0T2 ), although the spectra here start a 
day earlier, providing additional constraints. 

Regarding unburned material, other SNe la show dis- 
tributions of unburned carbon spanning a range of ve- 
locities. This suggests that carbon in some SNe la arises 
from incomplete burning, in contrast to the predictions 
of some delayed detonation models and merger scenarios 



( Marion et al.|2009| [Fryer et al.|201~0]|Ropke et al.|2012[ ). 

In general, due to the rarity of early SIN la observations, 
little is certain about the conditions under which carbon 
is observable in all subtypes. Noise makes weak carbon 
lines difficult to detect, especially when supernovae are 
faint at early times, and different velocities cause it to 
be blended with other features. Furthermore, SN la sub- 



types have different temperatures (Nugent et al. 1995), 
and these change with time, affecting whether carbon 
signatures exists as C I, C II, or C III. Thus carbon may 
be fairly common in SNe l a, but only apparent if spec- 
tra are obtained early (jParrent et al.|2011||Thomas et al. 
2011aHFolatelli et al.|20121|Silverman fe Filippenko|2012| 



Blond in et al.||2012p In addition, asymmetry may play 
a part in whether carbon is detected given that both 
explosive burning and WD mergers are inherently asym- 
metric processes (Livio & Pringle 2011 Pakmor et al 



2012 1. Spectropolarimetry and late-time observations, 
which are only possible for bright, nearby SNe like SN 
201 lfe, can shed light on these links. 

For the carbon in SN 201 lfe, we find C II to be 
strongest in the earliest spectra, suggesting it is mostly 
concentrated in the outer layers. The C II A6578 detec- 
tion persists until just before maximum light, although 
at later times it is very weak and would have been missed 
in any but the highest S/N spectra. Note, however, its 
presence throughout regions of freshly synthesized mate- 
rial shows that it is indeed distributed within the ejecta 
and is not just concentrated within the outermost layers. 
Because the high velocity O I is located above the outer- 
most reac hes of the C II, it is likely unburned progenitor 
material ( Nugent et al.j[20TT ). 

For the purposes of investigating spectroscopic diver- 
sity, it is convenient that SN 201 lfe exhibits the char- 
acteristic signatures of a "normal" SN la. Future com- 
parative studies of SNe la will benefit from these results 
and that of other well observed objects. The map of ele- 
ments in SN 201 lfe can now be used as input for model 
atmosphere codes, such as PHOENIX (see 



Hauschildt & 



Baron||1999 and references therein), enabling us to bet 
ter determine abundances, test theories of the explosion 
physics, and perhaps gain further insight into progenitor 
scenarios. 
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Figure 1. Optical spectra of SN 201 lfe during the first month of observations, labeled with respect to an explosion date of 2011 August 
UT 23.69 and a B-band maximum light date of 2011 September UT 10.3. Most of the data span 3500 - 9700 A, roughly every 2 days 
between the pre-maximum and post-maximum phases. Line identifications arc labeled for some of the most conspicuous features. The ©- 
symbol indicates the location of telluric features near 6279, 6880, and 7620 A. 
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Figure 2. Synthetic spectra comparisons are shown for days 1.5, 10.4, 20.6 and 32.6 post-explosion (—16.1, —7.2, +3.0 and +14.9 days 
with respect to maximum light). The dashed line represents the fitted continuum level. We color various pieces of our fits to indicate 
the dominant ion we attribute to the features. The black dotted line of the fit indicates where we do not attribute features to any of the 
assumed ions. 
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Figure 3. Symbols arc for clarity within the plot and different ions are represented by different colors, as labeled in panel (b). (a) Based 
on our time-series of synthetic SYNAPPS spectra, v m j n estimates of each ion are plotted versus time. The dashed line represents the location 
of our best modeled photosphere, (b) Here we reproduce panel (a), but include the Doppler widths of our fits. Color intensities are mapped 
to the normalized sum of a single absorption line component for each ion in every fit. Therefore, panel (b) only depicts where light has been 
scattered by a particular ion. Due to telluric features, we have only included our fitted Doppler widths for the high velocity O I component. 
The numbered rings (dashed- white) are in increments of 5 kkm s — . 
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Figure 4. "Snap-shots" of unburned material, i.e., spectroscopic signatures of unburned material as best traced by C II A6578 and possibly 
O I A7774. On the left, we compare our SYNAPPS fits (red line) to the observations (black line) that show a feature near 6300 A. At 1.2 days 
after t exp i the C II absorption minimum corresponds to a Doppler velocity of 16,000 km s _1 , at the upper reaches of photosphcric velocities. 
On the right, we plot the same but instead for the O I A7774 absorption feature. Our fits (blue line) place O I in both photospheric and 
higher velocity regions. The gray band highlights the early evolution of the O I feature and the ©- symbol indicates the location of telluric 
features. 



